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ABSTRACT 


A  FORTRAN  IV  program  which  implements  the  Polynomial 
Manipulation  System  (PMS)  is  presented  and  described.  PMS 
uses  the  Euclidean  Algorithm  to  reduce  a  system  of  poly¬ 
nomials  In  several  variables  to  a  resultant  system  which 
can  be  solved  sequentially  as  polynomials  in  one  variable 
(Kronecker's  method).  PMS  is  described  briefly  and  re¬ 
ferences  are  given  to  more  complete  discussions  and  to 
other  pertinent  literature. 
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I.  INTRODUCTION 


The  Polynomial  Manipulation  System  (PMS)  uses  the 
Euclidean  Algorithm  for  finding  the  eliminant  and  the 
greatest  common  divisor  (g.c.d.)  of  two  multi-variable 
polynomials.  All  polynomials  involved  are  represented 
symbolically;  PMS  is  a  computer  program  whose  Input  is 
the  symbolic  representation  of  two  polynomials  and  whose 
output  is  (normally)  the  symbolic  representation  of 
their  g.c.d.  and  eliminant. 

The  underlying  theory  and  the  application  of  PMS  to 
the  problem  of  solving  systems  of  polynomial  equations  is 
discussed  in  [^l"),  {[pj,  and[]^.  The  present  report  de¬ 
scribes  a  FORTRAN  IV  implementation  of  PMS  developed  on 
the  IBM  360  Model  50  at  Auburn  University. 

The  program  is  described  in  Section  II,  the  basic  flow 
charts  are  given  in  Section  III,  Input/Output  is  discussed 
in  Section  IV,  and  efficiency  of  the  method  is  discussed  in 
Section  V  along  with  possible  future  work.  The  FORTRAN 
program  is  reproduced  in  Appendix  A.  Appendix  B  contains 
a  simple  example  of  the  use  of  PMS  for  reduction  of  three 
polynomial  equations  in  three  variables  to  a  resultant  sys¬ 
tem  which  can  be  solved  in  sequence  as  polynomials  in  one 
variable . 
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II.  PROGRAM  DESCRIPTION 


The  PMS  program  is  basically  a  main  program  with  four 
subroutines,  only  one  of  which  is  significant.  The  other 
three  subroutines  are  used  for  output,  format  headings 
on  printed  output  and  scaling  of  coefficients  when  they 
become  large  enough  to  possibly  cause  an  overflow.  In  its 
present  form  the  program  is  limited  in  that  It  is  set  up 
to  use  only  175K  of  IBM  3^0  storage.  This  limitation 
places  constraints  on  the  program  which  allows  storage  of 
only  50,000  polynomial  entries  (each  term  has  n  +  1  entries 
where  n  is  the  number  of  variables  in  the  polynomial)  which 
are  presently  set  up  as  follows: 

1)  The  pair  of  polynomials  has,  at  most,  four  variables. 

2)  Each  polynomial  has  at  most  3160  terms. 

3)  The  leading  coefficient  polynomials,  to  be  defined 
below,  can  have,  at  most,  400  terms. 

Minor  modifications  could  Increase  the  number  of  terms 
or  variables  or  size  of  leading  coefficients  at  the  cost 
of  decreasing  the  others  or  by  use  of  a  greater  amount  of 
machine  storage.  Still  larger  polynomials  could  be  processed 
by  use  of  tape,  disc  or  other  storage,  but  this  has  not  been 
effected  since  such  increases  would  only  tend  to  accentuate 
certain  disadvantages  of  the  method  to  be  discussed  in 
Section  V. 

N 

Consider  the  pair  of  polynomials  U,T  :  E  ■♦R.  Let  x  . 

denote  the  variables.  Each  of  these  polynomial  functions 
can  be  considered  as  a  polynomial  in  x  whose  coefficients 


would  then  be  polynomial  functions  from  to  R.  Let  U° 

and  T°  denote  the  polynomials  In  x  which  are  the 

2  n 

leading  coefficients  of  U  and  T  respectively  considered  as 
polynomials  In  x^,  and  let  u  and  t  denote  the  degrees  of  U 
and  T  In  x^.  We  may  assume  t^u.  Consider  the  polynomial 
R  defined  by 

R  s  uOt  -  . 

1 

R  Is  a  polynomial  In  x^,...,x^.  Considering  R  as  a  poly¬ 
nomial  In  x^  with  polynomial  coefficients, It  Is  seen  that 
Degree(R)<t.  Let  Degree(R)  =  r.  If  r^u,  let  T  •»  R  and  re¬ 
peat  above  procedure.  If  r<u,  let  T  *  U  and  U  «  R  and 
repeat  the  above  procedure.  After  a  finite  number  of  ap¬ 
plications  of  this  algorithm  a  polynomial  R  will  be  found 
whose  degree  in  x^  is  zero.  Thus  R  will  be  a  polynomial  in 

X  ,x  ,...,x  .  It  is  easily  seen  that  at  each  stage  R  has 
2  3  n 

any  zeros  that  are  common  to  U  and  T.  The  R  which  is  free 

of  X  is  called  the  elimlnant  of  U  and  T. 

1 

III.  BASIC  PLOW  CHARTS 

The  flowcharts  for  output  and  scaling  will  be  omitted 
as  their  detail  Is  not  significant  to  the  main  purpose 
of  the  program.  The  main  program  flow  chart  is  given  on 
page  4. 
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MAIN  PROGRAM 


Read  U,T 

1 

Compute  U°,  T° 


Print  U,T 

i 

o 

I 

Determine  which  of  U  and 
T  has  greatest  degree  in 
Xj.  If  it  is  T  continue, 
if  not  interchange  U  and 
T,  U°  and  T°. 


Call  RESIDUE 
and  Form  R 


Scale  if  necessary 


\1^ 

If  R  is  free  of  x, ,  print 
R  and  return  to  beginning 
to  read  two  new  polynomials. 
If  not,  let  T  =  R. 


RESIDUE  SUBROUTINE 
Form  U^T 

Form  -T°Ux^“'^ 

i 

R  B  u°T-T°Ux^“'^ 

i  ‘ 

T  =  R 

I 

Calculate  new  T° 

I 

Is  t  =  0? 

/  \ 

no  N 

l/ 

Return  to 
Continue 
Processing 


yes 


Return, 
We  Have 
Eliminant 
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IV.  INPUT/OUTPUT 


The  polynomials  U  and  T  are  read  in  separately.  The 
first  card  for  each  polynomial  contains  the  number  of  vari¬ 
ables  in  the  polynomial  in  columns  31-3^  in  integer  format ^ 
right  Justified.  The  number  of  terms  in  the  polynomial 
appear  in  columns  35-38  in  integer  format,  right  Justified. 
Following  this  card  the  terms  of  the  polynomial  appear,  one 
term  per  card.  The  coefficient  appears  in  columns  1-16  in 
E  format;  following  this  are  the  exponents  of  the  variables 
right  Justified  in  integer  format  in  columns  17-21,  22-26, 
27-31,  etc. 

The  output  of  this  program  is  available  in  any  medium, 
altnough  the  program  is  currently  set  up  for  printed  output 
only. 


V.  EFFICIENCY  AND  FUTURE  WORK 

The  PMS  program  has  not  proved  useful  as  a  method  of 
reducing  polynomial  systems  of  equations  to  a  resultant 
system  for  the  following  reasons: 

1)  Storage  efficiency  is  low.  An  inordinate  amount 
of  core  is  needed  to  process  many  simple  appearing  systems 
of  equations. 

2)  Time  efficle..cy  is  low.  Extreme  amounts  of  time 
are  needed  to  solve  all  but  the  most  simple  problems.  As 
few  as  four  equations  in  four  variables  with  small  exponents 
(on  the  order  of  ten  or  less)  take  many  hours  of  machine 
time  to  reach  a  solution.  Simpler  problems  are  solvable  in 
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small  amounts  of  time,  but  other  methods  without  these  dis¬ 
advantages  can  be  used  to  solve  these  systems. 


3)  Certain  types  of  systems  give  solutions  which  have 
a  low  order  of  accuracy.  Several  articles,  PJ.  CsJ.  El. 
have  been  published  discussing  this  problem  as  well  as  the 
two  above. 

The  basic  PMS  program  will  be  examined  and  modified  to 
determine  if  it  is  of  value  in  algebraically  solving  simple 
systems  of  differential  equations. 
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AVPKW)TX  A 

FORTBA^^  IV  Prorrem  Listing 


r  )Xllk }((:({ 1^  ft  I9C  ft  jV y.  ft j'- ft:*::*-.  A:^  A  ft  ft  ft  ft  ft  ft  ft  A  ft  *  ft  ft  ):  ft  ftft  ft  ft  ft:  ?••.  :.-*:;;:  ft:  ft 

(  ftfti  ft.A  ft  ;•  A  ft  ft  ft  ft  ft  ft  ftft  ft  ?‘'ft  fxfi  5';  ft  ft  ft  s’-:;  A  ft  ;*:ft  •*•  :*t  ft  ?  /:  s';  .•':  :\i'z  •’  ft  ft:*  :'•  ft  ft  r’-.  ft  ft  ft 


I^'Pl  K, !  T  I-  '  :  =■  •'(  '  -•') 

f'V'.ir'.;  I  2  ./.•.'spi  ,  T  K-  (Arn  ) 

'U^,  0) ,  1 7'iA-:  ,  i  •  I  '  ,  (  M  -..1)  .  ^  ...}  ,  i 


V  !  'U  ?  1  f  0 ) ,  1 7  I  \ 

^ ;  s?,Lsi  ,1 

-  40  r 

(  *<.  ■-:  .■: 

I  •’  ;  •  :'r  ft  *  ft  ft  ft  ft  ft  <:  ft  ft  s',  ;  ft  r‘:  \ 


;  :*;  s’:  ft.  ft  ;;•  ft  .ft  : :  ft  ft  ::  <:  ft  s*  ft  ft  ft. :;:  ft  y.  ?••  ft  ;)•  V;  .'  ■  ft  ft  ft  ;?! 
ft-  ;•  *• ; :  -ft  ft.  ft  f/.  >\  i';  f  t  fi  ;  •  ft  ft  ;;-  ;•.  ft  ft  ft 


T  1  .. ,  ',  !  .  = 

I, '1  =  0 

»  7  t  "•• = 

W  A  (  ,1  ) 

»F(  'VArl. 

I  1  r.» 


P".!  Irl  ■  ■  I 

'  ■  i  (  I )  .  M  ■  (  !  .  .1)  ,  1  -  1 ,  OV  Ai:  i  • 
...(o,  I ' 

•’^■0.1  =  1,  t  •  : 

a;  .'  •)( 1  (  n  ,  (  1  ■  ( T ,  ,n  ,  1  :--i  ,nvak/) 

on  7  .1  =  It.iT  '  f- 1 

f’  f  .1  ?  I  ?  j ; 

(;M7[=i  ,  'iwa,  I 
5  n(  I  s  -  M 

IKNV/M'l  .Of.  -VAf’.?)  (  I  TO  106 

■lAX  =  WAS./ 

'J'/AHl  -■  1 

Oi’  ,!  -  1  ,.r,>; 

-,n  )  0  ’  :■  )  ,  -  I  -  ’:i 

-  i. 

f  \\/  I-,;;  1  n  ■ 

MAX  =  iO'Ai;  1 

jn^'ui  "  0 

ijn  'jis  ,!  =1,  ;T.  0V1 

IF(  .li’f.r,')  -  Iii(  )  ,.))  1  6-',-30,‘>0 

JflA;||  ^  HM  1  ,J  ) 

CO  .1  1  r;i  ip 
JU  ■■  ■ 

nn  j  -  i-.T'tf.-.i 

-  LM  1  ,  JS  )  05,34,33 

•  lU  =  .III  +  1 

KM' I'M  I!)  =  ;)(j) 

IF  I  A'.A'  I!0A  X(  .III  )  )  .1',..  1.F6)  LM  =  1 
or  ■-,■•1  ■-  :t  a,  max 


!j3  IUMAX(K-lf  JU>  =  IR(KfJ> 

55  -CONTINUE  -  - 

DO  10  J  =  l.NTERMl 

- - 

DOlOI=l,MAX 
10^  IU(  I  ,JT*rRUf  J) 

DO  9  J  =  1,NTERM2. 

RI JI=T< JI 
D09I=lfNVAR? 

9-  IRI  I,J)=I  Td  ,J) 

CALL  PR  i:'^T(2,NTERM2,0,0,NVAR2) 

- IRiNVAR2  .GE.  NVARl)  GO  TO  1075 

MAX  =  NVARl 
NVAR2  =  'JVAR2  +1 
on  104  J  =  NVAR2  iMAX 
00  104  K  =  1 ,NTERM2 
104  IR(J,K)  =  0 

NVAR2  =  MAX 
1075  MAX  =  NVAR? 

JPWRT  =0 
on  60  J  =ltMTPkM2 
IFtJPWRT  -  IR|l,Jn  51,60,60 
51  JPWRT  *  IK(1,JJ 

60  CONTINUE 

JT  =  n 

00  56  J  =  1,  NTfcRM2 
IFIJPWRT  -  IR(1,J))  56,57,56 
57  JT  =  JT  +  1 

TMAX(JT)  =  R(J) 

IF(AfiS( TMAX( JT ) I .GE. I.E6)  LM=1 
nn  50  K  =  2, MAX 
5»  ITMAXIK-l , JTI  =  IRIK.J) 

56  CONTINUE 

00  «  J  =  l,NTeRM2 

T( J)=R(J  ) 

0081=1, max 

R  IT( I  ,J)=IR{!,J)  , 

107  -  IFI  JPWRT-JPWRli)  70,71,71 

70  NN  =  JPWRT 

JPWRT  =  JPKRU 
JPWRU  =  NN 
MAXT=NTERMl 

IF{NTERM2.GT.NTFRM1)  MAXT=NTERM2 
00601=1, MAX 
TFMP=U( I ) 
ll<n=T{!) 

T(I)=TEMP 
OOeOJ=l,MAXT 
TEMP=IU( J, I ) 

!U(J,I)=IT(J,n 
flO  IT(J,I)=TEHP 

NN  =  JU 
JU  =  JT 
JT  =  NN 
NN  =  NTEPMl 
NTERMl  =  NTERM2 
NTFPM2  =  NN 
NN  =  NVARl 
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NVARl  =  NVAR2 
NVAft2  =  MN 
JJ=Jll 

TF< JT.GT.JJ) JJ=JT 
MA  =  MAX-l 
DO  73  J  =  1,  JJ 
TW=UMAX( J) 
t)MAX(J)  =  TMAX(J) 

TMAX( JJ=TW 
0073K=1,MA 
ITT=IUMAX(K, J> 

IIJMAX(K,J  )  =  ITf^AXlK.J) 

73  ITMAXlKtJ  n 

71  CONTIN'je 

C 

IP(lM,EO.l»  CALL  SCAlE2<NTrRMl,NTeRM?,MflX) 

LH  =  o 

101  CAll  PCSIOU  (MTFRMl,!vJTt«M2,MAX,J‘*ivRU,  JPWRT,  JU,  JT) 
IFdM  .F.O.  U  CALL  SC4LF2(MTtR'‘'l,NTHRM2,HAX) 
d(LM.^y.l)  G'Urt  3001 

IF(LS2.F(J.l)  CAIL  SCALE2<NTFRMl,NTl-k>^2,"AX) 

LS2*0 

3001  LM=0 

300  IFILSI.E'J.O)  GOTO  107 

L  S 1  =  0 

I  y=o 

301  CALL  I'RINT  (  3 ,  JU ,  NTFRH 1  ,  NTFR  .3?  ,  >•  AX ) 

HPITtCh,  mono)  ITIOCR 
10000  ,  OHSCALE,  13) 

GO  TO  100 
1  FORMAT! 30X , 21  A) 

4  rnRMAT(ElA.T,lOl3) 

FNO 


10 


C  SUPROviriNf  FOR  PRINTING  THf  TWO  POLYNOMIALS  AND  THE  RESIOJE 

C 

C 

SUeRHUTINF  PRINT  (L, JOfNTFRMl.NTFRMP.N) 

^^<1  4i  <:<!  #<1  <iii:  4^  4i4t<i  <i<t<t^]tc  4’^<<4>  4<4>4<4t4i4<4'  <'*« 

C 

IMPLICIT  INrEGFf<*2(  I-N) 

COMMIN  II  'AXI  400  I  ,  I  UMAX!  3 ,400  »  ,  1  T  MAX  (  3 , 400)  ,  TH  AX  (  400  )  ,  R I  BlfeO)  , 

NI  K(4,  Olf)''  )  ,  I  TPATL,  I  TlMFS,0{  31  60)  ,  IU(4, 31f  0)  ,T<  31C0)  -  IT(4  ,  3160)  , 

NL  S?,L  SI  ,L  1  I 

C 

f**  1,<«  <!  (;  IT  (t*.  4t 4141  •(;*[♦  1(. 4!  <141  4i 4< 41 1’/ 4t * 4i  4  4i4:  4> 4<* *  4i  4> 41  4> 4>4!  ti4t4<4t^4!4» 

C4>4<*4i*4<<*'('  +  «<<4v;/<'';'*i.''<'4'>'-<i*4'4>4<44<4'<'44>'<>fr4n(>4'*4>i;<4  ^♦X'  *4'4t  4>*^4'4'4>4<*4'4i4i4<4'4'4'4'<'4'4' 

IF(L.fO.-\)  TF  16, 310) 

I  f  (  I  .  '  C  .  P  )  ,^  '  I  T  c  (  6 , 3  1  I  ) 

ir(  L  ,F0.  1  )  WiU  TO  (  6,3)  .^) 


SO  TAIL  r-  rHlO) 

O'  I  3  I  6k  =  1  ,  J  ( I 

.'3IT:-t6,63  )  3{t.),nR{l,K),l-l,vj) 


31S 

CON  T  I  NO- 
H  =  N'-l 

503 

C  ■’N  1  I '  !U  ■ 

!\f  TO  •-'! 

43 

F0PMA1  (  1X,F  1^  .7,>y.,10n5,‘^X)  > 

311 

FnP*''AT(  IHI,' THH  C'OI  YNOMI  Al  0 

ISM 

31P 

FORM  Ml  1-1-,  '  TOP  f»i'l  YNOf*!  At  f 

IS*) 

310 

POR'.l  .\T(  )  ;-,  'THF  Ell^lMANI  IS* 

) 

540 

format  r-  jX  ,P  14  ) 

r>4i 

r  IPMAT  (  F  1  6  .  7,  1015) 

F\'0 

11 


c 

r 

r 

(, 

(, 

r, 

<’l 


?4 

n 


•iO 


■  SJB'^O'iTINE  FOR  PRINTING  hEAOINGS 
SUBROUTINE  PRir^THCK) 

THIS  SUBROUTINE  MERELY  PRINTS  COLUMN  HEADINGS  FOR  THE 

'  !  y<‘.riatles  depending  on  the  number  of  variables. 

■^HaT  is  its  only  purpose.  it  will  handle  up  to  10 
variables. 

GO  to  (  21, 22, 23, 2^,25, 26, 27, 2R, 29, 30)  ,•< 

WRITE(6,31) 

RETURN 

'  ■  WRiT,E(6<32) 

RET^/RN 

white  (6, '33)  i 
RETURN 

WHITE(6,34)  ! 

RETURN 

WRiTE(6,3S) 

RETI'RN^ 

WRITE|6,3U,)  ,  , 

RETURN 
WP.  ITEt6,37) 

RETURN  ■  '  . 

v.Hl  iE(6,  .-R  » 

RETURN 

WRITE(6,39) 

RETURN 
Wk  ITn  6,40) 

RET' 'N 

31  FORMAT!  IHO,  UHCOCF  f  IC  I E  NT  ,  I  OX  ,4HX  (  1)  ) 

32  FOR ''ATI  IHO,  UHCOEFE  IC  IF  NT  ,  lOX  ,  4HX  (  I )  ,  SX,4HX(2)  ) 

33  FORMAT!  lH0,llHCnEFFICIFNT,l0X,4HX(l),5X,4HX|2) ,5X,4HX(3)) 

34  FORMATI  1H0,11HC0EFFICIENT,10X,4HX(1),5X,4HX(2),SX,4HX(3),5X,4HX(4) 
1) 

35  FORMAT!  IHO , I IHCOEFF IC lENT , lOX ,4HX I  1) , 5 X, 4HX ! 2 ) , 5X , 4HX ! 3 ) , 5X ,4HX ! 4 ) 
1, 5X,4HX!  5) ) 

36  FORMAT! IHO, 1 IHCOF FF !C IFNT , 1  OX ,4HX 1 1 ) , 5 X, 4HX ! 2 ) , 5 X , 4HX ! 3 ) , 5X ,4HX (  4 ) 
1,5X,4HX(5) ,5X,4HX!6I ) 

37  FORMAT!  IHO , I IHCOEFF IC lENT , I  OX ,4HX ! I ) , 5X, 4HX ! 2 ) , 5X , 4HX ( 3) , 5X , 4HX I  4 ) 
l,5X,4HX!5) ,5X,4HX(6) ,5X,4HX!7)  ) 

38  FORMAT!  IHO , 1 1 HCOEFF I C IE  NT , lOX ,4HX ( 1 ) , 5X , 4HX ! 2 ) , 5X , 4HX ( 3) , 5X , 4HX ( 4) 
1,5X,4HXI5) ,5X,4HX(6),5X,4HX!7),5X,4HX|8) ) 

39  FORMAT! IHO, 1 IHCOEF F 1C lENT , lOX ,4HX ( I ) , 5X, 4HX ! 2 ) , 5 X, 4HX ! 3) , 5X ,4HX I  4 ) 
1,5X,4HX! 5) ,5X,4HX!6),5X,4HX!7) ,5X,4HX!8) ,5X,4HX!9) ) 

40  FORMAT!  IHO , IIHCOEFF I C If NT , I  OX ,4HX ! 1 ) , 5X, 4HK ( 2 ) , 5 X , 4HX ! 3) , 5X ,4HX ( 4 ) 
1, 5X,4HX! 5) , 5X,4HX(6) ,5X,4HX! 7) ,5X,4HX( 0) , 5X ,4HX( 9) ,5X, 5HX(  10) ) 

END 
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PEsr^ij':  stuRPUTiNP 


Sf'PTi,  TPMMS 


SURr^nUT  IN”  RFSI  Lnj(  NTU  tNTT  t^V  ♦  JPUf  JPT,  JU,  JT) 

C  *♦**♦*  <!*«**#**  <:  ♦<<#**************  1*1 !(.  ♦<!**  +  <"*+«>*%*<<♦*  -I'** 

Q**********  A*  **i,'i«ii(ii(<)!<**Di*<i*****i|i*<t«***i+****'^*i>«*^i>##*<<'l't**^  <'♦♦♦***<>*• 

c 

IMPLICIT  INT;GEf?*?n-N» 

(•.^M‘^V■'^l  I'"  AX(  40n  )  ,  rilMAXI  ^  ,400  )  ,  1  TM4V  (  3 , 400  )  ,  TM  AK  (  400  I  ,  R  (  3160  )  , 
MI  P(  4,  T1601  ,  1  T  vATP  ,  ni  Mt  S,IM  :3160  )  ,  IU(4 ,31  60)  ,T(  3lf>0)  ,  IT  (4,3160  ) 
NLSPtLSl.LM 

I)  IMF  NS  I  in  J:H5) 

MAXP=^lhO 

■lAXn  =  4  on 
TNTf'-.FR  Sk^ITCH 


C 

t  =1 

OO  2.<}  j  i,Lin 
A=T( J) 

1)04- Ktl,MV 

5210  1P(KI=1T(K,J) 

IP(  li  (  1 )  .OF  .,'PT  IcpTn^O 
|F(Ai>(A|  -  1."''') 

46  LS?=l 


RFTUfV'i 
•"'0  lO  L 
win  = 
IP!  i,n 


^  1,  .11) 

l'•■A.X(l)  >9' A 

=  f'ln 


lP<KF,n  -  H'MAX(nK-l,L)  ♦  l’.(KK) 
3M=l-l 

IF( MM.F  ).o  IG uTPlo 
nn80KK=l , MM 

IF(IRn,KK).NE.lR(l,!  nooToao 
r)!37nnKKK=2,MV 

I  F(  IKI  XKK,KK  J  .ME  .  IRIKKK,  I  )  )''.OTOPO 
CaNTINUE 

KCKX  )=q(k  K  )«-R(  I  ) 

1F(  APSCUKKJ  ).GT.1.F-9)GOT0  70? 

OT'OOILK-KK,NNM 

PI  IL)<)=P.  (  lLK+1  ) 

nN80  0NP;.)=  1  ,Ny 

IRl  mpo,  II  K  )  =  IPINPC,!  LK+1,  ) 

I  =1-1 
I  =  T-  1 
G3T "701 
('.ONf  I  MU' 

in  I  .LT.  'AX'’ )G0T()19 
WRIT:;(6,  12) 

FOPMATI  Dl  MANY  TEEMS  IN  Pi 

n  np 

!  =  1  +  1 
CONTI  'rj‘- 
nn  4.)  J  =  1,  JTii 

A  =  l)l.l) 


1  ’.(  KK  ) 


1.F-9)GOT070? 


KES  I  .JUE  ) 


I  B  <  K  I  -  I '  I  (  K  t  J  ) 

IK  Irt(  1  (.OF.  JPU  >001040 

1F(A^S(*>  -  47»4«t^Q 

4<^  LS?=1 

RFTOP.N 

4  7  Of)  41  L  =  It  JT 

f>(  I  )  =  -T-IAX  (|.)»A 
IP  (  1  , 1  >  =  OM  I  1  +  JPT  -  JPO 
Of!  3  0  ‘<  = 

30  IRIK, I)  =  IT3AX{K-1,L)  +  !«('<) 

MM=I-1 

)0  ^  bOKK  =  1 1  iV’ 

1  r  {  !•>(  I  ,KK>  .  'll;  .  Ik  (  1 ,  I  )  )  00103  F(' 
on^7onKKK=3,NV 

I  1-  (  1  :-•(  KI'F  ,KK  )  .NR  .  IP  (  KKK  ,  n  >001  (!'.p.0 

37C0  r.OijTlNO" 

■  (  K  K  )  S!K  ( '^  K  I  ♦  3(  I  ) 

I  (  A  -iS  (  '  (  K  1  )  .  (i  T  .  1  .P  -  9)  GCJTf)  3  f(:  (L 

M  0011  K=KK 

■V ( It  >=^' <  n K  +  i ) 


■)r'3'’  OO'I'-'I)-  1  ,  \'V 

3  M  Oi )  1 1'  (  ■■  t  1 1.  K  )  =  I  k  (  N  0 1  It  y  •*■  1  I 

1  =  1-1 
^707  ’=1-1 

fV^lO3  701 
7';!  f's-iIMUk 

1;  (  I  .t.T.  -’AX!  )G!:Tn41 
O'-  -11 

4]  ;  ^  I  +  1 


40 

IWI  !'■ 

T  r-  I  -  1 

99 

JO:,;;  T=0 

)  ^  -  0 

109 

1  -O^  ‘I  =  1,1 

I  F  (  I  9,  (  ]  ,  ’  )  -  JPHH-TI 

1 0  5,109  , 

10  ’ 

.|iO-l::T  = 

in-; 

00.9  I  !'\U- 
)  P  T  =  J  i’  'w  k  T 

I  ■'  (  '  P  .’F  T  )  1  3"  ,  1  a;  ,  1  ’.0 

1^-' 

1  SI  1 

lO  r-  I 

■  •;  7 y; 

1 

JT  =  0 
'■''’■•1  1  =  1  ,400 

T^■  A  <  (  1  I  )  =0  . 

O' ■  /  ,)  )  =  1  ,  T 

7 

IT'"  x(  j  1, 1  n  =0 

7F 

‘•(■'•'■•■F  J=l,  1 

1 :  (  JP/JRT  -  1  9.1  1  ,  J  )  ) 

9  0,  97,5.9 

)(  -  JT  1 

T^.-'KJTI  =  (J» 

1  r  (  T  'JI.AX  <  J  T  1 . 01-  .  1  .  f  <1 ) 
■rjc.)  K  =  .’,''JV 

L'-’=l 

1  T‘’.A  X(  K  -  1  ,  J  7  )  -  I  F  (  K 

,J) 

r 

V  riAjo- 

'  ;  l^.T  lAlll) 

n  Jl 

I  9  :• 


A 


DniJ=i,  I 
Tt J>sR( J) 
nniL=i,Nv 

1  ITIL,J»=1R<L,J) 

376  NTT=t 

RFTUR’') 


StlPRniJTINf;  SCALF2(  NTFRMl,;vlTeRM2,NV) 

IMPLICIT  •MTHGEI>*2(  I-NI  . 

COMMHM  U^AXI 400 ) , lUMAXI 3i400) f IT MAX ( 3,4001 ,TMAX(400 ) ,R( 3160) , 
NIR( 4, 11 ftO) , ITRATF, I T1 H6S , UI 3 160 ) , 1 U( 4 , 31 60) , T<  3160 ) , 1 T < 4 , 3 160 ) 
NLS2,LS1,IM 

FQIlTVAL6^!CE  (NVAR1,MAX) 

MAX  =  MV 

MVAK?sNVAR1  .  . 

on  I  I  =  1,400 

iJMAX(t)  =  U^AXt  I  ) /lOOO.  .  . 

1  TMAXd)  =  T^AXI  I  ) /lOOO. 

no  7  j  =  i,mtfkmi 
7  UIJ)  =  1J<J)/1000. 

on  9  J  =  i,ntcrm2  - -  - 

0  T(J)  =  T(.j)/1000. 

ITIMFS  =  ITIMFS  +1  - 

KFTURM 

FMD 
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APPENDIX  B 


This  appendix  presents  an  exainple  problem.  The  fol¬ 
lowing  system  of  three  polynomials  in  three  variables  is 
reduced  to  a  single  pol^Tnomial  In  one  variable: 

(1)  X  -+  X  +  X  =0 

12  3 

(2)  X  +  x2  =0 

1  2 

(3)  X  +  x^  =  0 

1  3 

First,  x^  is  eliminated  between  (1)  and  (2)  producing 
the  following  printout  which  has  been  labeled  for  expository 
convenience : 


COEFFICIENT 

O.IOOOOOOE 

01 

x(l) 

x(2) 

0 

x(3) 

0 

(1) 

O.IOOOOOOE 

01 

0 

1 

0 

O.IOOOOOOE 

01 

0 

0 

1 

COEFFICIENT 

O.IOOOOOOE 

01 

X(l) 

X(2) 

0 

X(3) 

0 

(2) 

O.IOOOOOOE 

01 

0 

2 

0 

COEFFICIENT 

O.IOOOOOOE 

01 

x(l) 

0 

X(2) 

2 

X(3) 

0 

(E  ) 

-O.IOOOOOOE 

01 

0 

1 

0 

-O.IOOOOOOE 

01 

0 

0 

1 

Second,  x 

is  eliminated 

1 

between  (3) 

euld  (2) 

producing 

the  following  printout: 

COEFFICIENT 

O.IOOOOOOE  01 

x(l) 

1 

x(2) 

0 

x(3) 

0 

(3) 

O.IOOOOOOE 

01 

0 

0 

2 

COEFFICIENT 

O.IOOOOOOE 

01 

X(1) 

1 

X(2) 

0 

x(3) 

0 

(2) 

O.IOOOOOOE 

01 

0 

2 

0 

COEFFICIENT 

O.IOOOOOOE 

01 

X(l) 

0 

X(2) 

2 

x(3) 

0 

(E,) 

-O.IOOOOOOE 

01 

0 

0 

2 
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Finally  (E^ )  and  (E^ )  are  treated  as  a  pair  of  poly¬ 
nomials  in  two  variables  x  and  x  .  Then  x  (x  in  our  first 

i  2  12 

system)  is  eliminated,  producing  the  following  printout: 


COEFFICIENT 

O.IOOOOOOE 

01 

X(1) 

2 

X(2) 

0 

(E,) 

-O.IOOOOOOE 

01 

1 

0 

-O.IOOOOOOE 

01 

0 

1 

1 

COEFFICIENT 

-O.IOOOOOOE 

01 

x(l) 

2 

x(2) 

0 

(E^) 

O.IOOOOOOE 

01 

0 

2 

COEFFICIENT 

O.IOOOOOOE 

01 

x(l) 

0 

X(2) 

4 

(Es) 

-0.2000000E 

01 

0 

3 

(Ep  is  our  eliminant  free  of  x^  and  x  ,  so  it  can 
be  solved.  Using  its  solution  (E^ )  can  then  be  solved. 
Using  this  (2)  can  be  solved  and  the  solutions  to  the 
resultant  system  (2),  (E^,  (E^)  are  the  solutions  to 
the  system  (1),  (2),  (3).  Three  other  equations  from  these 
six  could  have  been  taken  to  form  the  resultant  system. 
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